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ABSTRACT 


The  objective  of  this  study  is  to  analyze  the  behavior  of 
a  typical  ship  structural  member  exposed  to  fire. 

Included  in  this  work  are  chapters  pertaining  to  the 
modelling  of  the  fire  process  and  the  structural  member  to  be 
analyzed.  After  determining  the  expected  temperature  exposure 
of  a  shipboard  structural  member,  a  heat  transfer  analysis  is 
carried  out  and  the  resulting  load  is  subsequently  used  in  the 
failure  anal ys i s . 


The  major  thrust  of  this  study  is  the  analysis  of  the 
failure  mechanism  of  the  collapse  of  an  I  beam.  Two  methods 
ot  analysis  developed  by  Kecman,  and  Hayduk  and  Llierzbicki, 
are  summarized  and  app I i ed  to  the  case  of  a  fire  load.  The 
reduced  moments  due  to  fire  exposure  are  determined  and 
finally  recommendations  are  offered  for  the  use  of  this 
information  in  ship  structural  design. 
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Chapter  1 


TNiTPOm  IfTTTnM 


Tradi t i onal 1 y ,  Marine  Engineers  have  handled  the  effects 
of  fire  on  shipboard  structures  through  the  use  of  damage 
control  procedures.  Ideally,  a  structure  should  be  designed 
to  withstand  anticipated  fire  loads.  The  objective  of  this 
work  is  to  apply  the  methods  of  post  critical  analysis  of 
thin-walled  structures  developed  by  Mechanical  Engineers  to 
analyzing  fire  loads  on  board  ocean  engineering  structures, 
specifically  ship  structures. 

A  ship  structure  is  a  highly  redundant  structure.  If  the 
load  bearing  capacity  of  one  member  is  diminished,  then  the 
surrounding  structural  members  carry  more  of  the  applied 
load.  Shipboard  fires  are  char ac ter i st i cal  1 y  highly  localized 
phenomenon  which  may  effect  only  a  small  percentage  of  the 
entire  structure.  By  examining  the  effect  of  fire  on  a  single 
structural  member,  one  will  gain  insight  into  the  behavior  of 
the  entire  structure.  Ultimately,  ship  structures  may  be 
designed  to  withstand  the  effects  of  typical  shipboard  fires. 


It  is  the  intent  of  this  study  to  examine  the  effect  of  a 


•fire  on  a  typical  shipboard  structure,  namely  the  I  beam.  The 
major  thrust  will  be  the  analysis  of  the  failure  mechanism  of 
the  beam.  Two  methods  will  be  introduced  which  analyze  the 
post  yield  behavior  of  an  I  beam.  Figure  1.1  shows 
graphically  how  a  fire  load  diminishes  the  mechanical 
properties  of  a  beam. 

The  material  yields  upon  reaching  the  maximum  moment  in 
Figure  1.1.  Consider  the  following  relation: 

Stiffness  =  Momen t/Curvature . 

Beyond  the  yield  point,  stiffness  is  less  than  or  equal  to 

zero.  Normally  a  structural  member  remains  in  the  elastic 

range,  and  the  moment  remains  less  than  M  ,  the  fully  plastic 

o 

yield  moment.  However,  through  exposure  to  fire,  the 
stiffness  will  decrease  so  rapidly  that  the  moment  will  fall 
below  a  safe  range. 

The  ultimate  objective  of  this  study  is  to  show  how  much 
the  strength  of  an  I  beam  is  diminished  when  exposed  to  fire. 
To  this  end,  the  structure  and  the  fire  process  will  be 
modelled;  and  the  effect  of  temperature  on  the  mechanical 
properties  of  a  selected  material  will  be  given.  After  these 
models  have  been  developed,  a  heat  transfer  analysis  will  be 
performed.  As  a  result  of  this  analysis,  i t  wi 1 1  be  shown 
that  bending  collapse  is  the  primary  mechanism  of  structural 
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■failure.  Next,  the  -failure  mechanism  analysis  methods  will  be 
■fully  developed.  Finally,  recommendations  will  be  offered  for 
applying  the  results  of  this  work  to  the  design  of  an  entire 
ship  structure. 


Chapter  2 


Mnr>Fi  nr  RTPiirrruBE 


This  chapter  will  develop  the  motivation  towards  choosing 
an  I  beam  as  a  typical  shipboard  structural  element. 
Additionally?  a  material  will  be  selected;  and  the  temperature 
dependent  properties  of  that  material  will  be  given. 


2-J—Eajik.gcmiJnd 

The  principal  aspects  o-f  ship  structural  design  most 
commonly  employed  by  Marine  Engineers  involves  the  application 
o-f  rules  established  by  Societies  such  as  the  American  Bureau 
of  Shipping.  These  rules  are  based  upon  the  primary  loads 
experienced  by  a  ship  during  its  lifetime.  These  loads 
include:  still  water  bending  moments  and  shearing  forces;  wave 
induced  moments  and  shearing  forces;  impulse  loads;  thermal 
loads;  and  dynamic  loads. 


These  loads  which  are  unique  to  ocean  going  vessels  are 


often  difficult  to  accurately  predict.  Therefore,  ship 
structural  design  is  largely  a  process  involving  guesses 
derived  from  intuition  or  empirical  data.  The  primary  stress, 
strain  and  displacement  behavior  of  the  ship  structure  are  due 
to  displacements  limited  to  those  parallel  to  the  unloaded 
midsurface  of  the  plate.  Thus,  traditionally  hull  strength  is 
determined  through  an  analysis  of  the  midships  section. 
Figure  2.1  illustrates  a  typical  midships  section. 


Today,  Ship  Structural  Engineers  have  adopted  a  model 
which  is  motivated  by  considering  the  loads  transmitted  to  the 
structure.  Loads  transmitted  to  a  stiffened  plate  (by 
slamming,  etc.)  are  in  turn  transmitted  to  the  stiffeners  and 
from  them  to  the  point  supports,  and  so  on.  Eventually  a  load 
transmitting  grillage  is  derived.  Next,  boundary  conditions 
will  be  used  to  define  if  the  spans  are  treated  as  a  series  of 
isolated  points  or  continuous  beams.  In  most  cases  the 
boundary  points  will  be  those  which  can  be  described  as  simple 
or  clamped  supports. 

Figure  2.2  shows  examples  of  typical  members  of  a  ship's 
structure.  For  the  purposes  of  this  study,  an  I  beam  will  be 
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analyzed  as  its  behavior  is  considered  represen  tat i ve  of  the 
open  beams  in  Figure  2.2.  Figure  2.3  gives  the  dimensions  of 
the  beam  which  hereafter  will  be  used  as  the  model  to  be 
anal yzed. 


The  first  moment  of  area  of  the  I  beam  will  be  required 
in  the  analysis. 


b2t 

/  zdA  =  t  (b  -  t-)  (a  -  t  )  +  - - 

'  r  r  w  4 


(2.1  ) 


A  =  Section  Area 


t  »  Thickness  of  the  flanges 
f 

t  =  Thickness  of  the  web 
w 

a  ■  Width  of  flanges 


b  =  Height  of  web 

2 

Assume  t  is  negligible.  Hence,  the  first  moment  of  area  may 
be  approximated  by: 


pf  zdA  =  abt^ 


(2.2) 


2— 3_£alac±loa_o£_lda±aE.j.al 


Several  materials  have  been  used  throughout 
history  starting  with  wood,  and  progressing 
steel.  Today  several  grades  of  steel  and  even 


sh i pbu i  1  di ng 
to  iron  and 
t i tan i um  are 
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used  in  the  shipbuilding  industry.  One  of  the  most  common 
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grades  of  steels  used  is  ABS  Grade  A.  Material  Scientists, 
however,  have  not  tested  ABS  Grade  A  in  order  to  determine  the 
temperature  dependent  material  properties.  Ordinary 
structural  steels  such  as  ASTM-A36,  have  been  widely  used  as  a 
test  specimen.  Table  2.1  compares  the  properties  of  the  two 
steel s. 


Table  2.1 
Steel  Properties 


ASTM-A36 


ABS  Grade -A 


Carbon,  max  V.  0.23 
Phosphorous,  max  V.  0.04 
Sulpher,  max  V.  0.05 
Yield  Stress,  ksi  36.00 


0.26 

0.04 

0.04 

34.00  (min) 


It  is  considered  a  conservative  assumption  that  an 
analysis  for  ASTM-A36  can  be  used  without  sacrificing 
accuracy,  (ref.  1). 


One  of  the  major  contributions  to  the  loss  of  strength  is 
the  temperature  dependence  of  the  mechanical  properties  of  the 
materials  under  consideration.  Several  researchers  have 
predicted  the  deterioration  of  steel  structural  elements  at 
elevated  temperatures.  A  search  of  the  literature  reveals 
good  agreement  among  those  whom  have  studied  the  temperature 
dependence  of  material  properties. 
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Harmathy  and  Stanzak  (ref.  2)  have  carried  out  several 
experiments  in  order  to  provide  to  Design  Engineers,  primarily 
to  those  concerned  with  fire  endurance  of  structural  elements, 
the  elevated  temperature  mechanical  properties  of  common 
commercial  steel  products.  The  temperature  dependent 
expression  for  Young's  Modulus  of  Elasticity,  E,  (psi)  for 
ASTM-A36  steel  is  given  here: 

E  =  30*106  -  9.3T2  (2.3) 

o 

T  =  Temperature  of  heat  exposure,  Fahrenheit 

Steel  loses  strength  as  it  is  heated.  Figure  2.4  illustrates 
the  temperature  dependent  tensile  strength  of  ASTM-A3<6  steel. 

o  o 

Notice  that  between  1000  and  1200  F  the  tensile  strength 
falls  below  the  usual  design  stress. 

For  the  purposes  of  this  study,  linear  dependence  of 
temperature  is  assumed.  The  following  equation  will  be  used 
in  this  analysis. 

Cq  =  -20. 5T  +  40,000.0  4) 


1? 


I, 

Te^p^rAtur.e'  RvrvATiotns 
op  4STH-A36  "Ste-e^ 


F|  SUP.C  2.  i 


Chapter  3 


MnnEi  of  ft pf  PPnrFgg 


A  survey  of  the  literature  reveals  several  methods  of 
char  ac  ter  i  z  i  ng  a  -fire.  Harmathy  offers  a  numerical  technique 
which  is  summarized  in  Section  3.1.  Most  Fire  Engineers  use  a 
time  temperature  curve  to  define  the  characteristics  of  a 
fire.  The  development  of  an  acceptable  time  temperature  curve 
is  given  in  Section  3.2. 


3-»J_ £hacarl.eg.L3±xc.s-o^-a-ShjLjiboflrri  E-i-Cg 

The  heat  exposure  on  a  shipboard  structure  resulting  from 
a  fire  is  a  function  of  several  identifiable  and  incidental 
variables,  (ref.  3).  The  three  primary  factors  contributing  to 
the  severity  of  fully  developed  fires  is  the  flow  rate  of  air 
into  a  burning  compartment,  (compartment  ventilation);  the 
total  fire  load;  and  the  thermal  char ac ter i st i cs  of  the 
compartment  boundaries. 
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Numerous  theoretical  studies  have  determined 

t  i  me-temper  atur  e  curves  of  typical  -fires  based  on  these 
•factors,  <re-f  3).  Even  with  simplifying  assumptions,  these 
•factors  are  only  calculable  with  a  computer.  Harmathy 
introduces  three  easily  calculable  -fire  severity  parameters 
which  can  be  assumed  to  be  accurate  predictions  o-f  the 
contributing  factors  to  a  fully  developed  fire.  These 
parameters  are  defined  as  follows: 

1.  The  duration  of  the  fully  developed  fire,  t; 

2.  The  overall  penetration  flux,  q,  that  is  the  heat  flux 

absorbed  by  the  compartment  boundaries,  averaged 

spatially  over  the  boundary  surfaces,  and; 

3.  The  average  temperature  of  the  compartment  gases, 

(average  "fire  temper ature " >  ,  T  ,  averaged  spatially  over 

9 

the  period  of  full  development. 

It  is  useful  at  this  time  to  introduce  the  three  periods 
of  freely  burning  compartment  fires.  See  Figure  3.1.  The 
focus  of  this  study  is  on  the  period  of  full  development, 
since  80-90  percent  of  the  fuel  energy  is  released  during  this 
period. 


A  compartment  fire  is  a  fire  which  is  enclosed  in  a  shell 
of  arbi trary  shape.  Typically  on  board  a  ship,  ventilation  is 
secured  to  the  compartment  with  the  fire,  and  free 
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communication  to  the  surrounding  spaces  is  not  allowed.  Fire 
load  is  a  numerical  expression  which  is  a  function  of  the 
available  fuel  in  the  compartment.  The  units  of  specific  fire 
load  is  pounds  of  wood  per  square  foot,  <ref.  1). 

To  compute  the  fire  load  of  a  particular  compartment,  the 
products  of  the  weights  in  pounds  of  the  material,  and  the 
heat  capacity  are  summed  for  all  the  combustible  materials  in 
the  compartment.  The  resulting  expression,  which  has  units  of 
BTIT  s  must  be  divided  by  3000  BTU/lb.  This  result  will  provide 
the  same  available  energy  as  the  equivalent  weight  of  wood. 

Equation  (3.1)  gives  the  equivalent  weight  of  wood  in 
pounds  divided  by  the  deck  area  of  the  compartment. 

F  =  E  w . H . /H  A_  (3.1) 

ii  w  T 
i 

2 

F  ■  Specific  fire  load,  Ib/ft 

Ul  =  Weight  of  individual  materials,  lb 
i 

H  =  Heat  capacity  of  individual  materials,  BTU/lb 
i 

H  —  Heat  capacity  of  wood,  BTU/lb 
w 

2 

A  »  Area  of  deck,  ft 
F 

Table  3.1  gives  typical  fire  load  values  for  various 
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shipboard  compartments,  (ref.  4). 


Table  3.1 

Specific  Fire  Loading  of  Shipboard  Spaces 

lbs. /square  feet 

Control  Spaces 

Whee 1 house/Char troom  1.5 

Fire  Control  Stations  1.5 

Escape  Routes 

Corr i dors  1  . 5 

Stairway  Enclosures  1.0 

Accommodation  Spaces 
Staterooms 

Fire  Resistant  Furnishings  3.0 

Combustible  Furnishings  5.0 

Publ i c  Spaces 

Fire  Resistant  Furnishings: 

Lounges,  Restaurants,  etc.  3.0 

Ferry  Oesse Is  1.5 

Combustible  Furnishings  5.0 

Sanitary  Spaces  0.4 

Service  Spaces 

Gal  1 eys  10.0 

Pantries  (no  heating  appliances)  4.0 

Food  Concession  (ferry  vessel s- 

no  combustible  stowage)  1.5 

Workshops  10.0 

Storerooms  -  Combustibles  10.0 

Storerooms  -  Cleaning  gear  only  3.0 

Ship's  Laundry  10.0 

Laundry  -  Private  Use  1.5 

Main  Machinery  and  Cargo  Spaces  10.0 

Auxiliary  Machinery  Spaces  5.0 


Tanks  and  Voi ds 


0.0 


It  has  been  determined  by  several  researchers  that  the 
■flow  rate  of  air  into  a  burning  compartment,  (ventilation)  is 
just  as  critical  to  the  severity  of  the  fire  as  the  fire 
load.  The  determi nat i on  of  compartment  ventilation  depends  on 
how  the  rate  of  entry  of  air  compares  with  a  critical  value. 
The  rate  of  entry  of  air  into  a  compartment  can  be  expressed 
as: 


u 

a 


C,A  h 
1  w  w 


(3.2) 


U  =  Flow  rate  of  air,  Ib/hr 
a 

2 

A  =  Area  of  window,  ft 
w 

h  =  Height  of  window,  ft 
w 

5/2 

C  =  240  Ib/ft  -hr 
1 


The  critical  value  is  defined  as: 

(Ua!cr  =  C2*FAF 


(3.3) 


2 

C  =28.1  Ib/ft  -hr 
2 


w—  Material  constant  which  depends  on  the  furnishings  of 
the  compartment.  For  typical  living  spaces, 


0.55  <  9  <  0.90 


The  air  -flow  -factor  is  expressed  as  follows: 


(u  ) 
a  cr 


(3.4) 


Thus,  a  fire  is  ventilation  controlled  if  £  <  1.0;  and  a 
fire  is  fuel  surface  controlled  if  £  >=  1.0. 


The  average  rate  of  burning  of  combustibles,  R,  has  been 
determined  during  the  period  of  full  development  based  upon 
empirical  data.  For  £  >=  1.0 

R  =  C3<J>FAf  (3.5) 


2 

C  =4.57  Ib/ft  -hr 
3 


For  l  <  1.0 


R  =  0.163U 

3l 


(3.6) 


The  mass  flow  rate  of  gaseous  products  out  of  the 
compar  tmen  t  is: 

U  =  R  +  u,  (3.7) 

g  a 
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It  foil ows  that 


Q  =  R(0.9326H1  +  0.0686h2) 


<3. 8) 


where  : 


Q  =  average  rate  of  heat  evolution  from  combustible 
mater i al s 

H  =~  5760  BTU/lb,  heat  of  combustion  of  char 
1 

H  =~  11480  BTU/lb,  heat  of  combustion  of  volatile 
2 

decomposition  products 

S  -  the  fraction  of  heat  released  by  the  combustion  of 
the  volatile  decomposition  products  inside  the 
compar tmen t 

b  is  calculated  based  upon  an  estimation  of  the  average 
f 1 ame  length,  1  : 

l  =  c4(4)FAf)1/3  <3. 9) 

1/3 

C  *  1  .  1  1  f  t 
4 

*  =  c5ua1/3  (3.10) 

1/3  1/3 

C  *  0.366  f t(h)  /( 1b) 

5 
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Then  , 


if  £  =  H  6=1  <3.11) 

c 

if  £  >  H  5  =  (H  /£)1/3  <3.12) 

c  c 

where  H  is  the  height  o f  the  compartment, 
c 

The  duration  of  a  -fully  developed  -fire,  r,  is  calculated 
as  foil ows : 

for  £  =  1  T  =  C  /<$>  <3.13) 

6 

2 

C  =  0.205  ft  <hr )/l b 
6 

for  £  <  l  T  =  c  /*£  <3.14) 

6 

The  other  two  parameters  may  be  obtained  numerically  from 
the  following  equation; 

1  A  *  . 

q  =  7  (0  -  u  c  (CT  -  T  )  -  0-4(T  4  -  T  *)]  <3.15) 

A  *  g  g  g  o  3  g  o 

where  ; 

A  =  Total  area  for  the  compartment 
t 

Q  =  Average  rate  of  heat  evolution  <BTU/hr) 

o 

C  =  Specific  heat  of  compartment  gases  (BTU/lb-  R) 

9 


T  =  Temperature  of  combustion  gases,  R 
9 

o 

T  =  Temperature  at  time  =  0,  R 
o 

-8  2  o  4 

a  =  Stef an-Bol tzman  constant,  0.171*10  BTU/ft  -hr-£  R3 
£  =  Empirical  factor,  1.05,  dimensionless 


T 

g 


(T  +  f-(2TK/TT)1/2)  V/4 

O  K 


<3 .16) 


where : 

7  =  Empirical  -factor,  0.9,  dimensionless 

2 

«  =  Average  thermal  diffusivity,  ft  /hr 

o 

k  =  Average  thermal  conductivity,  BTU/ft-hr-  R 


Figures  3.2  and  3.3  show  the  dependence  of  the  overall 
penetration  flux,  and  average  temperature  of  gases, 
respectively,  on  the  air  flow  factor.  These  figures  were 

2 

determined  empirically  for  a  specific  surface  of  0.63  ft  /lb. 
Curves  are  given  for  three  values  of  the  specific  fire  load: 

2 

3.1 ,  6.1 ,  and  12.3  lb/ft  . 


Figure  3.2  shows  clearly  that  the  maximum  penetration 
flux,  and  fully  developed  fire  duration  occurs  at  [t  -  1  .  The 


specific  fire  loading  appears  to  have  little  effect  on  the 
maximum  flux,  and  the  duration  is  usually  quite  short.  It  may 
also  be  concluded  from  Figure  3.3  that  the  hottest  fires  occur 
in  ventilation  controlled  compartments.  It  is  possible  at 
this  time  to  draw  some  conclusions  about  shipboard  fires. 
Below  decks  fires  which  do  not  communicate  freely  with  the 
atmosphere  are  ventilation  controlled.  These  fires  are  more 
destructive  as  the  duration  is  greater,  the  average 
temperature  is  greater,  and  the  effective  heat  flux  is  greater 
than  for  an  above  decks  fire. 

McCarthy  demonstrated  the  usefulness  of  Harmathy's 
technique  for  predicting  fire  severity  of  a  typical  berthing 
compartment.  Figures  3.2  and  3.3  can  be  used  for  a  typical 
berthing  compartment,  and  in  accordance  with  Table  3.1,  a 

2 

specific  fire  load  of  3.0  lb/ft  was  used.  He  assumed  that 
the  ventilation  is  secured,  and  that  one  water  tight  door  is 
open  for  access  to  f i ref i gh ters. 

2 

A  =  15.0  ft 
w 

h  =  10  ft 
w 

U  =  11384.2  Ib/hr 
a 

2 

A  =  150  ft 
F 

®  =  0.75 
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Hence,  the  -fire  is  -fuel  surFace  controlled.  Notice  that 

as  A  increase  the  -fire  approaches  ventilation  controlled. 
F 

From  Figures  3.2,  and  3.3: 

3 

r  =  1.15*10  sec  =  19.2  min 

2 

q  =  6344  BTU/hr-Ft 
o 

T  -  1256  F 
9 

It  should  be  appreciated  -from  this  discussion  that  Fire 
is  indeed  a  complex  process  with  several  unquant i F i abl e 
variables.  The  term  destructive  potential  was  introduced  in 
reFerence  1.  The  destructive  potential  oF  a  Fire  may  be 
considered  a  Function  oF  the  maximum  temperature  obtained,  the 
duration  oF  the  Fire,  the  duration  oF  time  at  maximum 
temperature,  and  the  quantity  oF  energy  transmitted  From  the 
Fire.  Most  Fire  Engineers  use  a  time  temperature  curve  as  a 
deFinition  oF  heat  exposure.  For  the  purpose  oF  this  study, 
the  use  oF  an  accepted  t ime-temperature  curve  will  be 
considered  an  adequate  expression  oF  the  Fire  process. 


3->2-D^&lQfiman±,Q£-A-Iimj&-Xamfiat:a±una-£miu£ 


Traditionally,  fire  tests  have  been  performed  based  upon 
prior  established  standard  time  temperature  curves.  The  most 
widely  used  fire  test  in  this  country  is  the  ASTM-E11?  which 
was  developed  by  the  American  National  Standards  Institute, 
<ref.  5).  This  method  was  established  as  a  result  of  the 
necessity  to  standardize  the  structural  performance  of 
building  construction. 

The  basic  procedure  of  any  fire  test  is  to  place  the 
element  to  be  tested  into  a  gas  furnace  in  which  the  standard 
time  temperature  curve  is  duplicated.  The  gas  temperature  is 
closely  monitored  through  the  use  of  thermocouples. 
Additionally,  deflection,  axial  thrust,  and  collapse  behavior 
are  measured.  The  Ell?  time  temperature  curve  initially  rises 

o 

rapidly  reaching  1000  F  at  5  minutes  into  the  test.  Beyond 
this  point,  the  temperature  is  increased  more  slowly  so  that 

o 

at  the  end  of  one  hour  the  furnace  temperature  is  1700  F.  The 

o 

maximum  test  temperature,  2300  F  is  reached  at  5  hours  and 
held  constant  thereafter. 

Another  test  method  commonly  used  by  Ship  Structural 
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Engineers  is  the  International  Convention  -for  the  Safety  of 
Life  at  Sea  (SOLAS)  A60  test.  This  test  is  used  for  testing 
non-struc tur al  bulkheads  from  berthing  and  office 
compartments.  The  specimen  observed  is  not  loaded  but 
strained  at  the  edges.  Other  non-struc tural  bulkhead  tests 
are  SOLAS  A30 ,  A15,  AO,  B15,  BO,  but  the  requirements  for 
these  tests  are  not  as  severe  as  for  the  A60  test,  <ref.  1). 

For  the  purposes  of  this  study,  reference  6  supports  the 

o 

use  of  a  step  change  to  1700  F  vice  the  gradually  increasing 
temperatures  of  ASTM-E119  or  the  SOLAS  A60  curves.  Figure  3.4 
gives  the  ASTM-E119  curve,  the  SOLAS  A60  curve,  and  the 
assumed  time  temperature  relation  to  be  used  in  this 
analysis. 
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Chapter  4 


1 


I 


The  intention  of  this  analysis  is  to  determine  the 

relation  between  temperature  and  the  distance  along  the  length 

of  the  beam.  A  beam  exposed  to  a  fire  load  will  be  modelled 

as  a  fin  with  a  temperature  gradient.  The  temperature  at  the 

exposed  face  is  T  ;  the  temperature  at  the  end  of  the  fin  is 

b 

T  ,  the  ambient  temperature.  The  geometry  of  the  fin  is 

a 

defined  in  Figure  4.1 


The  following  assumptions  are  used: 

-  T  =  Temperature  at  x  —  0 

b 

-  T  =  Temperature  at  x  =  L 

a 

-  The  heat  transfer  is  steady  state. 
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-  k  =  Uniform  conductivity  of  the  material, 

o 

BTU/C in-sec-  F] 

-  h  =  Uniform  surface  coefficient  of  heat  transfer, 

2  o 

BTU/C hr-ft  -  F3 . 

Additionally,  this  study  will  use  the  assumption  that  the 
rod  is  thin,  which  implies  that  the  temperature  in  the  thin 
direction,  perpendicular  to  x,  is  uniform.  Therefore,  the 
surface  temperature  equals  the  centerline  temperature,  (ref. 
7).  Applying  an  energy  balance  on  an  elemental  section  *x , 
yields  a  governing  differential  equation.  In  steady  state, 
with  no  internal  heat  generation: 


Sc  Sc+Ax  +  ^Aa 


(4.1) 


Expanding  the  first  term  of  the  right  hand 
following  equation: 


Sc  ~  Sc+Ax  =  Sc  '  (Sc 

-  \ 

Sc  ~  Sc+Ax  dT*” 


da 

+  ■^_~Ax) 


si  de 


gives  the 


(4.2) 

(4.3) 


Sc  Sc+Ax  ■^•(-kAyAz'l^-)  Ax 
Sc  Sc+Ax  4—  (ks4~)Ax 


(4.4) 


(4.5) 


Where  S  =  Surface  area  =  AyAz  .  Conduction  heat  transfer  across 
A A  is  given  in  the  following 
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equation: 


a 


q.  =  h(T  -  T  )A 
^Aa  a 

q.  =  hPAx (T  -  T  ) 
Aa  a 


Where 


P  =  Perimeter  =  2<Ay  +  Jz) 

T  =  Temperature  of  fin  at  x 

Substituting  <4. 5)  and  <4. 8)  into  <4.1)  yields: 

^(kS^)Ax  -  hPAx  (T  -  TM  =0 

Finally,  the  governing  differential  equation  for  0 

T  ]  is: 
a 

Where  S,  k,  and  T  are 

a 

The  differential 
state  heat  flow  with  one  internal  heat  generation  source 

A  - 

dx2  k 


—  -£|e-o 

,  2  kS 
dx 


constant. 


equation  for  one  dimensional 


<  4 .  <£) 

(4.7) 

<4. 8) 


<4. 9) 

=  [T  - 

<4.10) 

steady 
i  s : 

(4.11) 


Note  that  the  fin  problem  is  similar  to  the  single  internal 


generation  source  problem  with  a  negative  heat  source: 


W 


i 


hP9 

S 


<4.12) 


4-.2_Solulj.on 


The  method  of  solution  will  -follow  a  method  provided  in 
reference  8.  Assume  the  solution  is  o-f  the  form  given  in 
equat ion  <4.13): 


9 


Cje8*  ♦ 


(4.13) 


Where  : 


B  =  [hP/kS]1/2 


<4.14) 


Tne  boundary  conditions  for  a  fin  of  finite  length,  L,  as 
defined  by  Figure  4.1  are: 

@  X  =  0;  9  =  Tb  "  Ta  5  eb  <4.15) 

@x=L;  6  =  T-  T=  0  <4. 16) 

a  a 

Substituting  <4.15)  and  <4.16)  into  <4.13)  yields: 


C1  +  C2 


BL  _  -BL 
Cle  +  C26 


<4.17) 

<4. IS) 
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Solving  for  C  : 

1 

a  <rBL 

c  =  "ebe  (4.1?) 

1  BL  -BL 
e  -e 

C  follows  from  (4.17).  Finally  the  solution  is: 

2 

e  eB(L-x)  +  e-B(L-x)  (4.20) 

6~  BL  -BL 

b  e  -  e 


At  this  point,  it  is  necessary  to  solve  for  B.  For  the  I 
beam  in  Figure  2.3,  the  perimeter  to  surface  area  ratio  is: 


Figure  4.2 
temperature 


P  2(b  +  2a  -  t  ) 
—  =  w 
S 


bt  -  2t ,t  +  2t,a 
w  f  w  f 


(4.21 ) 


gives  the  relation  between  thermal  conductivity  and 
for  plain  carbon  steel,  (ref.  7).  Within  the  range 
o  o 


of  interest,  70  F  <=  T  <=  2000  F,  the  curve  is  considered 

approximately  constant.  For  the  purposes  of  this  analysis, 

o 

the  thermal  conductivity  at  1000  F  wi 1 1  be  assumed. 


k  =  0.0005  BTU/in  sec  °F 


(4.22) 


The  condition  of  natural  convection  in  a  gas  is  assumed 
for  this  analysis,  (ref.  7): 

h  =  2.0  BTU/hr  ft2  °F  (4.23) 

Therefore , 

1/2 

B  =  (P/10.8S)  '  (4.24) 

where,  P/S  is  given  in  equation  (4.21) 
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The  validity  of  thin  fin 
According  to  reference  7,  the 


theory  must  now  be  checked, 
thin  fin  theory  is  valid  if: 


hb  J 


(4.25) 


Substituting  equations  (4.22)  and  (4.23)  into  (4.25),  the 
maximum  value  for  b  is: 


b  =  3.6  ft 


(4.26) 


This  is  a  reasonable  criteria  for  this  study.  Therefore  thin 
fin  theory  is  valid. 


4*3_Sasui±s 

It  is  now  possible  to  draw  some  conclusions  about  the 
heat  transfer  through  the  beam.  A  Microsoft  BASIC  computer 
program  was  written  in  order  to  examine  several  geometric 
configurations  of  the  beam.  The  source  code  for  the  program 
is  in  Appendix  A.  The  results  are  given  graphically  in  Figure 
4.3.  The  dimensions  of  the  beams  examined  are  in  Figure  4.4. 

The  following  conclusion  may  be  drawn  from  Figure  4.3.  As 
P/S,  and  L  increase,  the  effect  of  the  fire  becomes  more 
localized.  With  a  typical  shipboard  I  beam,  the  ratio  of 
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perimeter  to  surface  area  is  such  that  the  fire  load  may  be 
treated  as  a  point  load.  That  is  to  say  that  the  heat 
transfer  through  an  I  beam  remains  highly  localized  in  the 
vicinity  of  the  fire.  All  effects  of  the  fire  virtually 
disappear  at  a  distance  which  is  a  quarter  of  the  length  of 
the  beam  away  from  the  fire. 


Chapter  5 


ANALYSIS. QF-IHF ■EAliiiEE  MECHANISM 


As  was  stated  in  the  preceding  chapter,  an  I  beam  exposed 
to  a  typical  shipboard  -fire  wilt  experience  the  same  -failure 
mechanisms  as  a  beam  under  a  point  load.  Recently  there  has 
been  increasing  interest  in  this  type  o-f  -failure  due  to 
mechanical  loads.  It  is  the  intention  o-f  this  study  to 
approximate  the  shipboard  -fire  load  as  an  additional 
mechanical  load,  and  to  use  two  o-f  the  recently  developed 
techniques  -for  analyzing  the  failure  mechanism.  Additionally, 
this  chapter  will  include  summaries  of  literature  which 
focuses  upon  bending  collapse  of  thin  walled  structures. 

Kecman  (ref.  10)  studied  the  bending  collapse  behavior  of 
rectangular  and  square  section  tubes  both  theore t i cal  1 y  and 
experimentally.  He  used  a  limit  analysis  technique  and  a  set 
of  formulae  to  derive  the  relation  between  the  hinge  moment 
and  the  associated  angle  of  rotation.  Uery  good  agreement  was 
found  between  his  theoretical  predictions  and  the  experimental 
results.  His  theoretical  predictions  were  in  the  form  of  a 
numerical  solution.  However,  Hayduk  and  Wierzbicki  (ref.  11), 
gave  a  closed  form  analytical  solution.  In  reference  11, 


lower  and  upper  bound  solutions  for  the  mean  crushing  strength 
of  cruc i forms  were  obtained  by  considering  bending  and 
extensional  modes  of  deformation. 


Kecman  developed  a  numerical  scheme  to  analyze  the 
bending  collapse  of  rectangular  and  square  section  tubes, 
<ref.  10).  The  primary  failure  mode  is  local  buckling  and 
hinge  development.  He  considered  only  uniaxial  bending 
col  1  apse . 


5.1.1  Maximum  Bending  Strength 

Kecman  referred  to  a  study  which  examined  the  bending 
strength  of  a  box  shaped  section  with  the  webs  in  tension  and 
the  top  flange  in  compression.  The  effect  of  the  webs  was 
found  to  be  negligible.  Therefore,  the  critical  stress  of  an 
I  beam  may  be  assumed  to  be  the  same  as  the  critical  stress  of 
a  box  shape  beam. 

The  critical  stress,  a  ,  of  the  section  was  derived  in  a 

cr 


5 


I 


■form: 


0  =  - - 2-( 5.23  +  0.16^)  £) 

Cr  12 (1-V  )  b  a 

a  =  -f  1  ange  wi  dth 


<5.1 ) 


b  =  web  height 

t  =  web,  and  -flange  thickness 

I-f  the  critical  stress  is  less  than  the  yield  stress, 

the  material,  then  the  -flange  buckles.  The  effective 
becomes : 


a  Of 
o 

width 


ag  =  a  (0. +  0.3) 
e 


<5. 2) 


o  -  stress  at  the  edge  of  the  flange 
e 

If  the  critical  stress  is  greater  than  or  equal  to  the  yield 

stress,  then  a  =  a.  For  the  purposes  of  this  study,  assume 

e 

that  collapsing  occurs  when  a  =  a  .  I t  follows  that  the 

e  o 

maximum  moment,  M  ,  may  be  calculated  using  <5.1)  and  <5. 2). 

m 


If  a  <=  a  : 
cr  o 


2a  +  b  +  a  (3-  +  2) 
e  n 


M  =  0  tb  * - 

m  o  3  (a  +  b) 


<5. 3a) 


If  »  >=  3  o-  : 

cr  o 


Mm  *  M0  =  °0t[a(b  -  t)  +  0.5(b  -  2 1 ) 2 ] 


<5. 3b) 
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where  M  =  the  -fully  plastic  yield  moment  o-f  the  section.  If 
o 

a  <  a  <  3  <j  ,  a  moment,  M  is  de-f  ined: 
o  cr  o  o 

M  1  =  o  tb  (a  +  — )  <5.3c> 

o  o  3 

With  the  resulting  maximum  moment! 

a  -  a  _ 

M  =  m  '  +  (m  -  m  1 )  r _ 2.>  ^ 5 . 3d) 

m  °  o  o  2a  1 

o 

5.1.2  Actual  Bending  Collapse  Mechanism 

A  detailed  analysis  o-f  the  hinge  deformation  pattern 
offers  insight  into  the  failure  mechanism.  The  deformed 
pattern  is  given  in  Figure  5.1. 

The  following  charac ter i st i cs  of  the  deformation 
mechanism  were  given  in  reference  10: 

<a)  Most  of  the  plastic  deformation  is  concentrated  along 
the  stationary  yield  lines:  JF ,  HI,  EB,  68,  FC,  HC,  CR, 
BA,  GK,  and  EL. 

<b>  Travelling  yield  lines  of  rolling  deformation  were 
also  observed  along  GA,  AE,  KA,  and  LA  while  the  distance 
BA  changes  during  hinge  rotation. 

<c)  The  corner  angle  along  KL  remains  almost  constant. 


<d)  The  deformation  of  the  tension  flange  includes  the 
travelling  yield  lines  NP,  and  MG,  appr ox i mate  1 y 
correspondi ng  to  the  intersection  of  the  undeformed 
flange  planes  and  the  current  cylindrical  surface  TU.  See 
Figures  5.1,  and  5.2. 

<e)  In-plane  deformation  of  the  walls  was  also  observed, 
particularly  along  the  rolling  lines  GA,  AE,  BA,  HC,  CF, 
IR,  and  JR.  This  deformation  is  essential  to  maintain 
hinge  kinematics. 

<f>  The  rolling  radius,  r,  varies  slightly  along  GA, 
measured  by  cutting  at  right  angles  to  GA  and  its 
estimated  value  was  approximately  proportional  to  the 
hinge  length. 

r  <  0 .03-0 . 05) *KL 

The  radius  reduces  during  coll  apse ;  so  that  rolling 
becomes  more  and  more  difficult,  until  the  mechanism  jams 
and  starts  a  new  phase  of  collapse. 

<g)  The  rolling  radius  along  KA,  and  LA  varies  from 
infinity  at  K  and  L  to  a  value  close  to  r  in  <f)  above. 

Four  distinct  phases  in  the  hinge  development  were 
den  t i f i ed : 

<1)  Initial  phase  with  a  bulge  in  the  web,  A  ,  in  Figure 


<2>  Second  phase  with  rolling  as  described  by  the 


cross-section  shape  A  in  Figure  5.3. 

2 

<3>  Third  phase  indicated  by  a  jamming  of  the  rolling  and 

described  by  the  cross  section  B  A  D  in  Figure5.3. 

3  3 

<  4)  Final  phase  initiated  by  the  contact  o-f  the  two 
buckled  halves  o-f  the  compression  flange.  There  is  no 
clear  transition  from  one  phase  to  another. 


5.1.3  Theoretical  Bending  Collapse  Mechanism 

Since  structures  most  commonly  experience  hinge  collapse 
described  by  the  second  phase  which  includes  rotation  angles 
in  the  range,  5-10  <=  <=  25-35.  Kecman  based  his  theoretical 

analysis  on  the  following  assumptions: 

-  The  web  and  the  flanges  of  the  section  deform  along  the 
concentrated  yield  lines  only; 

-  The  web  and  flanges  are  incompressible  and 
inextensible; 

-  The  structural  continuity  is  maintained  in  two 
char ac ter i st i c  sections;  and 


-  Some  in  plane  deformation  is  essential  for  hinge 
k i nemat i cs. 

The  model  is  given  in  Figure  5.3. 

The  rolling  radius  is  negligible  in  comparison  with  the 
other  hinge  dimensions.  Points  A,  5,  C,  and  R  are  in  a 
straight  line.  The  coordinates  of  point  B  are: 

B(x,y,z)  =  (h,  bcosp  -  [bsinp(2h-  bsinp)]1^2,  0)  (5.4) 

From  continuity  of  the  middle  section,  the  web  height,  b,  is 
de  term i ned : 

,  2  2.1/2 

b  -  zft  +  <yA  +  *A)  (5.5) 

With  y  -  y  ,  and  equations  (5.4),  and  (5.5): 

A  B 

z  =  bsin2p  -  hsinp  +  [bsinp(2h  -  bsinp)  1  *cosP  (5.4) 
A 

An  analysis  of  the  rotation  of  the  fiber  in  the  top 

flange  will  lead  to  the  equation  for  the  hinge  length,  2h .  For 

a  fiber  which  is  initially  parallel  to  the  beam  and  passes 

through  point  A,  longitudinal  continuity  is  implied.  The 

distance  from  this  fiber  and  the  edge  KL  i s  AD  =  b  -  z  . 

A 

Initially,  the  coordinates  of  the  fiber  and  the  line  GK  are 

CO,  b-z  ,  03.  After  rotation,  the  coordinates  are  C(b 

A 

z  )sin/>,  (b-z  )cos  p,  03.  Point  A  has  the  coordinates  Ch,  y  , 
A  A  A 
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I 

I 


s 


z  ].  The  distance  from  point  A  and  the  rotated  point  is  h  if 
A 

the  web  and  flanges  are  i nex tens i bl e .  The  equation  for  h  is 
theref ore  : 

h2  *  [h- <b-zft) sinp] 2  +  [yA~ (b-zA)cosp] 2  +  z  ^  (5.7) 

Substitution  of  <5. 4),  and  <5. 6)  into  <5. 7);  and  meeting  the 
condition  required  by  <5. 5),  it  may  be  said  for  any  ratio  a/b, 
and  any  rotation  angle: 

2h  =  a  or  2h  =  b  <5. 8) 


The  value  which  yields  the  lowest  energy  absorption  will  be 
used.  Thus, 

< 

2h  =  a  if  a  =  b 


2h  *  b 


if  b  *  a 


(5.9) 


The  second  phase  of  deformation  is  complete  when  the 

rotation  angle  equals  the  jamming  angle,  a  .  Jamming  starts 

J 

when  the  two  buckled  halves  of  the  compression  flange  contact 
each  other.  From  Figure  5.3  the  jamming  angle  is: 


a  =  2arcsinp  (-  — - 

J  D 


<5. ICO 


5.1.4  Nominal  Energy  Absorbed  During  Hinge  Rotation 


The  assumptions  listed  in  the  previous  subsection  imply 


that  all  the  energy  is  plastic  and  is  absorbed  along  the  yield 
lines  of  the  theoretical  collapse  mechanisms.  Relative 
rotation  along  the  flange  edges  JF,  and  IH  results  from  hinge 
rotation,  see  Figure  5.1.  From  geometry,  of  Figure  5.3: 

7T  t 

Y  =  j  -  P  -  arcsin  (1  -  j^sinp)  (  5 . 1 1  > 

The  energy  absorbed  along  these  edges  becomes: 

W1  =  WJF  +  WIH  =  2mpa[^-  -  p  -  arcsin (1  -  jj^sinp)  ]  <5.12) 

m  s  Average  fully  plastic  moment  of  a  unit  width  of 
P 

the  flange 


Similarly,  it  may  be  said: 


w„ 


=  V.’  =  m  a  {  tt  - 

CR  p 


2arcsin(l  -  r-sinp)  ] 
n 


<5.13) 


The  energy  absorbed  by  the  yield  line  AB  can  be  found 

from  <5.11),  where  the  length  of  AB  =  z  is  known  from  <5. 6), 

A 

and  the  relative  rotation  is  <  *  -  2/5).  Hence: 


W  =  w  =  m  [bsin2p-hsinp+ (bsinp  (2h-bsinP) )  1//2*cosP]  [iT-2arcsin  (1-^sinD)  ] 
3  AB  p  n 


Similarly,  along  BG,  the  length  is  h,  and  the  relative 
rotation  is  n/2 : 


w„  =  V*  +  VJ  =  m  hn 

4  BG  BE  p 


<5.15) 


Figure  5.4  shows  the  relative  rotation,  *1  ,  along  yield 


K.  c? 


The  Anklet  of  Rotation 


FicwiR*  5,4 


1 


line  GK.  The  construction  of  this  figure  is  detailed  as 

follows:  the  plane  co  is  normal  to  GK  and  contains  point  A. 

xy 

Point  A'  is  the  projection  of  A  onto  the  xy  plane.  A  is 

o 

obtained  by  rotating  o>  about  its  intersection  w  with  the  xy 

xy 

plane.  Coordinates  of  A"  are: 


1/2 

htanp  +  bcosp  -  [bsinp (2h-bsinp) ] 


V  = 


1  +  tan* 


XA"  =  yA,,tanP 


<5.16) 


<5.17) 


The  energy  absorbed  along  both  the  GK,  and  EL  1 i nes  are  then: 


Wr  =  w_„  +  Vtf  =  2m  barctan  [ - — - - 

5  GK  EL  p  l(h-xs..)2  +  (ys..-YB)211/2 


(5.18) 


The  rolling  radius  along  the  yield  lines  GA,  and  AE 
remains  relatively  constant,  so: 


W_  =  t»  +  W„„  =  2m  (-)z 
6  GA  AE  p  r  A 


<5. 1?) 


During  hinge  rotation,  the  rolling  radius  varies  as 


described  by: 


r  =  (0.07  -  ^)h 


<5.20) 


Equation  <5.20)  is  based  on  empirical  data.  Kecman 
approx imated  the  variati on  of  the  rolling  radius  al ong  KA  with 
a  hyperbolic  approximation: 


'KA  £ 
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<5.21  ) 


Where  1  is  the  distance  from  K  along  KA.  The  linear  variation 
k 

of  the  rolled  length  between  K  and  A  is  given  by: 


o  =  J£z 

r  KA  A 


<5.22) 


The  energy  along  KA  is: 


2m  z  KA 
P  A 


<5.23) 


W=W  4.  W  _  «  Ar,  ^  l  Z. 

7  WKA  +  WLA  "  3  p~7th  +  ^  +  ZA  1 


2.1/2 


<5.24) 


The  relative  rotation  along  yield  lines  PN,  and  QM  i s  p  , 
whereas  along  KL,  and  MN  the  angle  is  given  by: 

ZA 

C  =  arctan  ( — ) 
yA 

Finally,  the  energy  absorbed  by  these  yield  lines  is 
approximated  by: 


w8  =  WPN+WQM+WKL+WMN  =  2mp[ap  +  Cretan  (—)  ] 

J  A 


<5.25) 


At  this  point,  the  total  energy  absorbed  by  the  hinge  may 
be  expressed  as: 


w(a)  =  I  w. (a) 

i=l 


<5.26) 


This  result  is  valid  for  a  rotational  angle  less  than  the 


jamming  angle.  Kecman  determines  the  moment-angle 


relationship  numerically  -from  the  -following  equation: 


M(a) 


w(a  +  Aa)  -  W(a) 
Aa 


(5.27) 


Where  represents  a  small  but  -finite  increment  o-f  the  hinge 
angle,  a. 


This  -formulation  assumes  that  the  material  is  ductile 
enough  not  to  -fracture  during  hinge  rotation  or  subsequent 
bending  collapse.  This  is  considered  a  safe  assumption  for 
common  structural  steels.  The  average  plastic  moment  is  given 
by  the  following  equation: 


m  = 
P 


(5.28) 


The  yield  stress,  and  the  maximum  nominal  flow  stress  as 
functions  of  temperature  were  determined  in  Chapter  5. 


5.1.5  The  Final  Form  of  the  Hinge  Moment 

The  analysis  in  the  preceding  section  is  not  valid  for 

the  first  phase  of  deformation.  This  phase  will  be 

approximated  by  assuming  that  for  a  moment  greater  than  the 

fully  plastic  moment,  M  ,  is  constant  and  equal  to  M  .  This 

o  o 

would  be  a  conservative  assumption.  For  hinge  angles  greater 
than  the  jamming  angle,  Kecman  offers  the  following  formula: 

M(a)  =  M(Oj)  +  1.4 [M  -  M(aJ)](a-  oij)  (5.2?) 
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This  -formulation  has  been  programmed  in  a  M  i  crosof  t-BASI  C 
program.  The  source  code  is  given  in  Appendix  B.  Figure  5.5 
shows  the  results  o-f  the  program. 


The  primary  di-f-ference  between  Kecman's  method  and  this 
method  is  that  Hayduk  and  Wierzbicki  considered  not  only 
i nextensi onal  deformation  modes,  but  also  extensional 
deformation  patterns.  The  effect  of  these  extensional  modes 
will  be  fully  developed  in  this  section. 

Hayduk's  and  W i erzb i ck i ' s  formulation  results  in  a  closed 
form  analytical  solution  which  predicts  the  mean  crushing 
force.  Reference  11  specifically  addresses  the  problem  of 
analyzing  the  crushing  process  of  thin  walled,  plate  -formed, 
open  structures  with  particular  emphasis  on  L  shapes  and 
cruci forms.  The  I  beam  can  be  approximated  by  two  T  shapes  or 
four  L  shapes.  Lower  and  upper  bound  solutions  for  the  mean 
crushing  strength  of  cruciforms  was  obtained  by  considering 
various  modes  of  deformation.  According  to  Hayduk  and 
Uierzbicki,  extensional  deformation  accounts  for  at  least  one 


third  of  the  dissipated  energy 


The  end  products  of  this  section  will  include  the 
application  of  Hayduk's  and  Ul  i  er  zb  i  ck  i  '  s  method  to  an  I  beam, 


and  reviews  of  both  references  11  and  12. 


5.2.1  Collapse  Modes 

Hayduk  and  Wierzbicki  considered  out  o-f  plane 

deformations  in  which  the  deflections  exceed  10-100  times  the 

section  thickness  and  are  of  the  order  of  magnitude  of  the  web 

height  and  the  flange  width.  The  process  of  collapse  starts 

as  the  compression  flange  loses  stability,  either  elastically 

or  plastically  at  a  critical  load,  P  .  Following  the  buckling 

c 

deformation,  the  flange  can  maintain  added  loads  until  a 

maximum  load,  P  ,  is  reached.  Eventually,  as  the  resistance 
max 

decreases,  the  structure  collapses  completely.  The  structures 
being  examined  here  often  regain  strength  after  the  initial 
loss  of  resistance.  Subsequent  load  peaks  reveal  the 
development  of  additional  folds. 


5. 2. 1.1  Isometric  Transformation  Modes  < I nex tens i onal 

Deformat i ons) 

These  modes  are  of  the  type  used  by  Kecman  to  analyze  the 
collapse  mechanism.  In  this  type  of  deformation,  the  energy 
is  concentrated  over  hinge  lines,  while  the  rest  of  the 
structure  remains  rigid.  An  assumption  of  i nex tens i onal 
deformation  is  that  material  continuity  is  maintained.  During 
the  transf ormat i on  of  two  surfaces,  the  lengths  of 
corr espond i ng  arc  elements  is  preserved.  Hayduk  and 
Ulierzbicki  show  the  progression  of  various  deformation  paths. 
These  are  illustrated  in  Figure  5.6. 

In  Figure  5.6  <a),  and  <b>,  the  final  state  is  formed 
i nex tens i bl y .  This  can  only  be  done  if  the  boundaries  have 
sufficient  freedom  to  deform.  The  element  in  Figure  5.6  (c) 
is  initially  folded  along  the  stationary  hinge  line  AB; 
finally  the  element  is  rotated  extensibly,  <ref.  11). 


5. 2. 1.2  Hinge  Collapse  Modes 

There  are  two  types  of  hinge  collapse  modes.  The  first 
is  a  propagating  hinge  line  which  occurs  when  the  hinge  line 
is  free  to  move  with  respect  to  material  points,  and  the 
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material  of  one  -flange  is  e-f  f  ec  t  i  ye  1  y  transmitted  to  another 
•flange.  Extensional  de-formation  occurs  in  the  vicinity  o-f  the 
hinge  line.  This  mode  is  a  char ac ter i st i c  of  thin  walled 
structures  consisting  of  two  intersecting  plates  such  as  the  I 
beam  being  considered  in  this  study. 

The  second  hinge  collapse  mode  is  the  formation  of  a 
stationary  hinge  line.  This  mode  is  char ac ter i zed  by  a  hinge 
line  fixed  in  the  material.  Side  flanges  which  are  initially 
rectangular  are  transformed  into  trapazoidal  elements. 
Therefore,  extension  occurs  over  an  entire  flange.  This  mode 
will  not  be  considered  in  this  analysis,  (ref.  11). 


5. 2. 1.3  Assembled  Collapse  Modes 

Hayduk  and  WierzbicKi  demonstrated  the  method  of  assembling 
the  folding  patterns  shown  in  Figure  5.6.  The  patterns  are 
folded  into  more  complicated  symmetric,  asymmetric  and  mixed 
collapse  modes  which  are  shown  in  Figure  5.7. 


I  Moving  Hinge  Un 


Pe TFORKATlOH  R\YT£-^ 
Figure"  5.  3 


5.2.2  Anal ys i  s 


For  the  purposes  of  this  study,  the  analysis  of  a  T 
section  will  be  used;  this  analysis  is  considered  to  be 
equivalent  to  an  I  beam  analysis.  In  reference  12,  the  short 
T  cross  section  under  consideration  was  subjected  to 
compression  with  clamped  boundary  conditions.  A  mode)  is 
presented  in  Figure  5.8. 

For  the  purpose  of  estimating  the  mean  crushing  force, 

P  ,  the  T  section  is  expanded  into  a  plate  and  an  angle  as 
m 

shown  in  Figure  5.9. 

The  plate  will  collapse  through  the  bending  zone  limited 
by  the  two  horizontal  moving  lines.  The  angle  will  follow  the 
folding  mechanisms  described  by  Figure  5.10(b).  The  collapse 
modes  in  Figure  5.10(b)  is  consistent  geome tr i cal  1 y  and 
kinematically.  The  regions  are  described  as  follows: 

(I)  Four  plane  trapazoidal  elements  as  rigid  bodies. 

(II)  Two  sections  of  cylindrical  surfaces  at  which 
continuous  bending  takes  place  without  any  extension. 

(III)  Two  sections  of  conical  surfaces  in  which  material 
is  bent  and  rebent  again  as  the  material  moves  from  one 
flange  to  another. 


Plate 


Artc,u= 


FoU>lHO,  [i 

-Fi^wrte  5. 10 


(IV)  A  section  of  toroidal  surface  which  produces 

extension  in  a  c i rcumf ren t i al  direction  and  continuously 
changes  principal  curvature  in  the  other  direction,  (ref. 
11) 


The  following  assumptions  will  be  used  in  this  analysis: 

-  Neglect  the  short  side  of  the  angle  element. 

-  The  length  of  the  two  horizontal  hinge  lines  is 
approx imatel y  equal  to  the  length  of  the  large  side  of 
the  angle  element. 

-  Two  inclined  hinge  lines  are  of  a  Tength  equal  to  a 
half  wave  length,  H. 

-  The  length  of  the  toroidal  surface  in  the  vertical 
direction  is  neglected. 

The  energy  dissipated  of  an  angle  element  is  given  by 
summation  of  three  terms,  (ref  11): 

1.  The  energy  of  a  continuous  deformation  field  in  the 
section  of  a  toroidal  shell: 

E1  “  9-28M0Hr/tw  (5.30) 

2.  The  energy  dissipated  by  the  discontinuous  velocity 
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■field  in  the  horizontal  hinge  line  o-f  length,  Is 


E_  =  l 
2  2  o 


<5.31) 


3.  The  work  done  by  an  inclined  hinge  line  through  the 
deformation  process: 

E3  =  2.22MQH2/r  <5.32) 


where : 

r  =  The  small  radius  of  the  toroidal  surface. 

2 

M  *  cr  [  t  ]  /4,  the  fully  plastic  yield  moment  per  unit 
o  o  w 

width 

H  =  Mode  half-wave  length 


a  —  0 . 7<t 
o  ou 

a  =  Ultimate  strength 
ou 

t  ■  Web  thickness 


w 


For  the  models  given  here: 


2 HP  =  E,  +  4  (E_)  .  +  4(E  +  2E 

m  1  2  web  2  flange  3 


<5.33) 


After  substitution  of  the  dissipated  energies  into  <5.33): 


t*  2 


VMo  =  4*64f  +  ^[c  +  d(~)  J  +  2.22^ 

w  «  w 


<5.34) 
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Optimizing  equation 

(5.34)  by  using  the  following 

partial 

der i vat i ves : 

1 

( 5 . 35a) 

3P m 

3H  “  0 

(5.35b) 

Equating  (5.35a)  and 

1  (5.35b)  yields: 

r  = 

,2.22  2  .  .fcf  2H1/3 

'  2  ™„ <c  *  d'r> 

4  •  64  w 

(5.36) 

H  = 

■n2t.  t  , 

(2.22)  (4.64)  (c  +  d(~>  )  1 

w 

(5.37) 

Ultimately,  the  three  basic  mechanisms 

of 

energy 

dissipation  appear 

equally  weighted  in  the  final 

equation  for 

mean  crushing  force: 

P  /M  = 
m  o 

t  2 

»  3  1  (4.64)  (2.22) tt(1  +  k(^)  )  (~)  ]  ' 

w  w 

(5.38) 

Equation  (5.38)  may 

be  simplified  as  follows: 

p  / 
m 

Mo  =  3 [AlA2A3 (|-)  ] 1/3 
w 

(5.39) 

where : 

k 

c 

(5.40) 

A!  =  4.64 

(5.41 ) 

t-  2 

a2  =  +  k(~)  J 

(5.42) 

w 

a3  =  2.22 

(5.43) 

Note  also  that  equations  (5.34),  (5.36),  and  (5. 

37) 

may  be 
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simplified  into  the  following  forms: 


■1  = 
o 

Al~  +  A2H  +  A37 

W 

(5.44) 

r  = 

[<A„A,/A.2)  (ct  2)]V3 
-3  X  W 

in 

^r 

in 

•w 

H  = 

t (a|/aiA  ) (c2t  ) ] 1/3 

(5.46) 

The  application  of  these  equations  to  the  case  of  an  I 
beam  requires  that  the  geometric  parameters  be  modified.  An  I 
beam  is  two  T  sections,  hence  the  dimensions  are  doubled;  but, 
b  =  2c,  which  yields  the  following  relations. 

/M  —  *3  r  7\  *  *  /  V,  /hi,  \  1  1/3 


P/M  =  3fA.A-A, (b/2t  )]' 
mo  12  3  w 


k  =  2a/b 


(5.47) 


(5.48) 


It  is  now  possible  to  apply  this  relation  determined  for 
the  crushing  case  to  the  bending  case.  The  geometry  of  this 
crushing  mode  is  given  in  Figures  5.11  and  5.12. 

Recall  the  propagating  hinge  line  mode  as  described  in 
paragraph  5. 2. 1.2.  This  mode  describes  the  global  geometry  of 
the  crushing  process  and  is  determined  purely  by  the  width  of 
the  two  side  panels,  d,  which  is  prescribed,  the  wavelength, 
H,  which  is  to  be  determined,  and  the  crushing  rotation  angle, 
n,  which  is  a  parameter  of  the  process.  The  relation  between 
the  axial  displacement,  f> ,  and  the  angle,  <*,  is: 


<S  =  2H  ( 1  -  cos3) 
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Reference  12  determined  an  expression  for  the  instantaneous 


crushing  force,  P,  normalized  with  respect  to  P  s 


p/p»  '  it2c°s6  +  SKe  +  sln‘B> 


<5.50) 


Figure  5.13  gives  the  geometry  of  the  bending  case.  From 


Figure  5.13: 


a  6/2 

tan2  *  b 


<5. 51) 


Thus,  it  may  be  said  that: 


P/P  =  f (6)  =  f(2btar£) 
in  2 


<5.52) 


For  a  section  width,  b: 


H-|b 


<5.53) 


M  =  —  b 
m  2 


<5.54) 


Where  M  is  the  mean  moment  due  to  bending  rotation,  a,  as 


defined  in  Figure  5.14.  After  equating  <5.53),  and  <5.54): 


M _ P 

M  _  P 
m  m 


<5.55) 


Figure  5.15  gives  P/P  as  a  function  of  <V2H  based  on 


equations  <5.49),  and  <5.50),  <ref.  12).  This  graph  is  an 


equivalent  description  of  M/M  as  a  function  of 


< b/H) £  tan<«/2) 3 . 


pEFiNi'Tiwrt  of  An^I£  of  Rotation 
Figure:  S.  l*i 


For  this  purposes  of  this  study,  an  approximate  relation 
■for  M/M  and  a  Mi  1  1  be  determined.  Assume  equation  <5.56)  is 


a  valid  -function  for  M/M  : 


M/M  =  D/(6/2H)n 
m 


<5.56) 


where : 


D,  n  are  constants 


Assume  the  new  curve  goes  through  points  P  ,  and  P  ,  in  Figure 

1  2 

5.15. 


P,  (<5/2H,  M/M  )  =  (0.30,  1.0) 
J-  m 


<5.57) 


P,  (<S/2H,  M/M  ) 
*  m 


=  (1..0,  0.24) 


<5.58) 


Solving  for  D,  n: 


D  =  0.24 


<5.59) 


n  =  1.2 


<5.60) 


Assuming  n  =  1.0,  the  following  equation  results: 


<5.61 ) 


Figure  5.16  shows  how  the  approximate  solution  compares  to  the 
exact  solution.  Substituting  equation  <5.51)  into  equation 


<5.61),  the  following  results: 
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<5.62) 


M/M  _  0124 _ 

7  m  -  (b/H) tan (a/2) 


Considering  equation  <5.54): 


(b/H) tan (a/2) 


<5.63) 


Now  substitute  equation  <5.47)  into  <5.63) 


» -  ^vva'^v)173’  tb'21 


<5.64) 


But  M  is  the  -fully  plastic  yield  moment  for  the  flange 
o 

a  t2 

v-V- 


<5.65) 


where  t  =  t  =  t  . 

f  w 


The  desired  final  result  of  this  analysis  is  a  relation 

between  M  and  M  ' .  Rewriting  equation  <5.65): 
o 

M  '  =  a  /  zdA 
o  o  A 

Recall  equations  <2. 2)  and  <2. 4): 

CMT)  =  -20. 5T  t-  40,000.0 

2 

J  zdA  =  abt  + 

a  4 

Finally,  substitute  equation  <2. 4),  <5.66),  <5.65),  and 

<5.46)  into  <5.64): 


<2. 2) 

<2. 4) 
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A  .7  -3  -H  r  .6  .7  -S 


Results  °f  H  ArfDuv^  ah\p 
^/eR2-6icKi$  Analysis 
F\G,U9.£  s.YI 


M/M  ’  = 
o 


3a  t 
o 


4aQ (abt  +  b  t/4) 


0.24 [(A^/A  A  )b22t]1/3 
[A1A2A3 (b/2t) 1  {2] [  btan (a/2)  5 


The  simplified  final  result  of  this  relation  is: 


M/M  ' 
O 


0.36TTt(b  +  2a) _ 

tan(a/2) [b(b  +  4a)] 


<5.68) 


Figure  5.17  shows  graphically  the  resulting  moment-angle 
rotation  relationship  as  determined  by  Hayduk  and  Wierzbicki. 
Appendix  C  contains  the  source  code  for  a  Microsoft  BASIC 
program  which  determines  the  moments  resulting  from  equation 
<5.72)  . 


The  primary  difference  between  the  two  methods  of 
analyses  is  that  Kecman  considered  only  i  nextensi onal 
deformation  modes,  whereas  Hayduk  and  Wierzbicki  considered 
extensional  deformations  as  well  as  i nex tens i onal 
deformation.  Since  the  second  method  considered  more  energy 
dissipation  than  the  first,  the  allowable  moments  due  to 
deformation  are  greater.  However,  for  angles  greater  than  one 
tenth  of  the  jamming  angle,  Kecman's  method  yields  greater 
moments.  The  errors  induced  by  his  numerical  integration 
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technique  are  compounded  as  the  angles  increase.  Figure  5.18 
shows  graphically  how  the  results  of  each  analysis  differs. 
The  dimensions  of  the  section  analyzed  are  given  below: 

a  =  12  inches 

b  =  1 2  i nches 

t  =  t  -  .6  inches 
f  w 

The  maximum  moment  given  in  Figure  5.18  is  the  fully  plastic 

yield  moment  of  the  section,  hereafter  called  M  ,  which  is 

o 

only  a  function  of  geometry  and  temperature  of  exposure  of  the 

beam.  Hence,  M  is  not  a  function  of  the  analysis  method, 
o 

Clearly,  as  the  angle  of  rotation  increases  the  moment 
decreases;  and  as  the  temperature  increases,  the  moment 
decreases.  Appendices  B  and  C  also  contain  results  which  show 
the  effects  of  changing  the  geometry  of  the  section.  For  both 
analyses,  the  moments  are  decreased  as  the  section  thickness 
decreases;  the  web  height  decreases;  or  the  flange  width 
decreases . 

These  methods  of  analyses  are  considered  accurate 
predictions  of  the  failure  mechanisms  of  a  beam  exposed  to 
fire.  Recommendations  for  using  the  information  provided  by 


these  analyses  will  be  given  in  the  following  chapter. 


Fisuae  5  l ft 


G 


Chapter  6 


r.ntjr  i  t  is  i  n *un  PEcnMMFNDAT  I  QMS 


A  summary  of  each  chapter  is  presented  here  with 
concluding  remarks.  Finally,  recommendations  are  offered  -for 
using  this  analysis  for  future  studies. 


■d-J—  Summary 

As  stated  in  the  introduction,  the  ultimate  objective  of 
this  study  was  to  show  how  much  the  strength  of  a  shipboard 
structural  member  is  diminished  when  exposed  to  fire.  To  this 
end,  Chapters  2,  3,  and  4  provided  input  into  the  analyses  of 
the  failure  mechanisms  studied  in  Chapter  5.  The  results  of 
the  preceding  chapters  will  be  summarized  here. 

Chapter  2  developed  the  model  of  the  structure.  Simple 
observations  revealed  that  the  ship  is  a  redundant  structure 
composed  of  beams,  plates,  and  columns.  The  motivation 
towards  choosing  an  I  beam  for  the  analysis  was  developed. 


The  materia]  selected  for  the  model  was  ASTM-A36  steel  as  its 
physical  and  mechanical  properties  are  similar  to  the 
properties  of  ABS  Grade  A  steel,  a  common  shipbuilding  steel. 
The  equations  for  the  temperature  dependent  mechanical 
properties  were  the  output  of  this  chapter  to  be  used  in  the 
analysis. 

Chapter  3  provided  background  for  char ac ter i z i ng 

shipboard  fires.  Three  parameters  which  are  assumed  to  be 
accurate  predictions  of  the  contributing  factors  of  a  fully 
developed  fire  were  defined.  These  factors  are:  the  duration 
of  the  fully  developed  fire;  the  overall  penetration  flux;  and 
the  average  fire  temperature.  Since  fire  is  indeed  a  complex 
process  with  several  unquant i f i abl e  variables,  a  time 
temperature  curve  was  used  as  a  definition  of  the  fire 
process.  For  the  purposes  of  this  study,  a  step  function  to 
o 

1700  F  was  used  as  the  time  temperature  relation. 

Chapter  4  was  the  heat  transfer  analysis  which  determined 
the  relationship  between  the  elevated  temperature  exposure  and 
the  distance  along  the  length  of  the  beam.  The  results  of 
this  analysis  led  to  the  conclusion  that  a  beam  exposed  to 
fire  will  experience  the  same  failure  mechanisms  as  a  beam 
under  a  point  load. 

Finally,  the  primary  thrust  of  this  study  was  to  analyze 
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the  -failure  mechanisms  due  to  fire  exposure.  Two  methods  were 
fully  developed  in  the  preceding  chapter  and  the  results  were 
summarized  graphically  in  Figure  5.13. 


The  ultimate  objective  of  any  structural  design  process 
is  to  develop  a  structure  which  is  strong  enough  to  withstand 
all  expected  service  loads  and  a  high  percentage  of  unexpected 
loads.  The  information  provided  by  the  analysis  in  Chapter  5 
is  considered  useful  input  for  t.<e  design  of  shipboard 
structures.  It  is  possible  through  the  use  of  Figure  5.18  to 
determine  how  the  stiffness  of  a  structural  member  is 
decreased  when  exposed  to  elevated  temperatures .  The 
diminished  stiffnesses  may  be  used  as  input  to  a  finite 
element  program  which  simulates  a  ship  structure.  The  effect 
of  a  shipboard  fire  upon  the  entire  structure  may  then  be 
de  termi ned. 

The  examination  of  the  effect  that  one  weakened  member 
has  on  the  other  members  of  a  redundant  structure  is 
considered  fertile  ground  for  future  studies.  Reference  16 
may  provide  insight  into  the  behavior  of  this  phenomenon. 
Zubelewicz  and  Mroz  numerically  simulated  rock  burst  processes 
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treated  as  problems  of  dynamic  instability.  According  to 
reference  16,  local  plastic  failure  occurs  at  one  location 
(such  as  in  a  member  of  a  ship's  structure)  and  the  elastic 
energy  flows  to  the  hinge  and  is  the  driving  force  into  the 
damage.  Energy  then  flows  into  the  neighboring  hinges,  which 
is  to  say  that  one  unstable  member  may  produce  a  chain 
reaction  and  other  members  may  subsequently  fail. 

As  a  final  word,  it  is  highly  recommended  to  utilize  the 
information  provided  in  this  study  as  an  aid  for  examining  the 
effect  of  fires  on  the  ship's  structure.  Ultimately  a  ship 
may  be  designed  to  withstand  shipboard  fires. 
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Appendix  A 


a-EEOGEAM 


10  REM  PROGRAM  TO  CALCULATE  TEMPERATURE  GRADIENT  IN  A  BEAM 
20  DIM  TEMP<20> ,Y<20) ,2(20) ,TOP<20) 

25  INPUT  "ENTER  THE  BEAM  LENGTH  IN  FEET" ; XL 
30  INPUT  “IS  BEAM  SOLID  <S>  OR  I  < I > " ;N* 

40  IF  N*="I"  THEN  GOTO  90 
45  LPRINT  “  SOLID  BEAM- 
50  INPUT  “ENTER  WIDTH  OF  BEAM" ;W 
60  INPUT  "ENTER  DEPTH  OF  BEAM" }D 
65  LPRINT  “  WIDTH,  DEPTH®" ;  W;D 
70  B=SQR<<12/5.4)*<<1/D)+<1/W>>> 

80  GOTO  155 
84  LPRINT  "  I  BEAM" 

90  INPUT  "ENTER  WIDTH  OF  FLANGE" ;W 
100  INPUT  "ENTER  DEPTH  OF  WEB" ;D 

105  LPRINT  "  WIDTH  OF  FLANGE,  DEPTH  OF  WEB*";  W;D 
110  INPUT  “ENTER  THICKNESS  OF  FLANGE" ;T1 
120  INPUT  "ENTER  THICKNESS  OF  WEB" ;T2 

125  LPRINT  "  THICKNESS  OF  FLANGE,  THICKNESS  OF  WEB®" ; T 1 ;T2 
130  P®2*<  D+2*W+2*T1 -T2) 

140  S=2»< T 1 #W) ♦ < D#T2) 

150  B®SQR<  < 1 2/5 . 4) »<  P/S) ) 

155  LPRINT  "  LENGTH  OF  BEAM®"; XL; "  FEET" 

156  LPRINT  "  B»";B  s  LPRINT  "  " 
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r 


. 


160 

BL=B*XL 

D 

170 

BOT=EXP <  BL  > -EXP  < -BL ) 

; . 

■ 

• 

180 

X—XL/20 

■ 

■ 

190 

LPRINT  “  PT  NO  X  T-Ta/Tb-Ta" 

200 

FOR  KX=1  TO  20 

210 

Y<KX)=X*KX 

220 

Z<KX)=B*<XL-Y<KX>> 

6 

230 

TOP<KX>=EXP<Z<KX) >-EXP<-Z<KX>> 

240 

TEMP  <  KX  >  =TOP  <  KX ) /BOT 

250 

LPRINT  "  " ;KX,Y<KX) ,TEMP<KX> 

m 

260 

NEXT  KX 

k 

270 

END 

B 


M 


i 


r 

\ 

t 

L 

[ 

r 

L 


— j 


i 

i 
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Appendix  B 


10  DIM  UK 20) 

20  INPUT  "TEMPERATURE  OF  EXPOSURE,  DEGREES  F“ ;T 
30  INPUT  “FLANGE  WIDTH,  INCHES" ;A 

40  INPUT  “WEB  HEIGHT,  INCHES" ;B 

41  INPUT  "THICKNESS,  INCHES" ;TH 

42  LPRINT  "  ■ 

50  LPRINT  "  " 

51  LPRINT  "  TEMPERATURE*" ;T 

52  LPRINT  “  FLANGE  WIDTH*"; A 

53  LPRINT  "  WEB  HEIGHT**'  ;B 

54  LPRINT  “  THICKNESS*" ;TH 
60  S0*-20 . 5*T+40000 . 0 

65  IF  T>=1 951  THEN  S0*4.5 

67  MO*  SO  *B*TH*  < A+  B/4 ) 

68  LPRINT  "  FULLY  PLASTIC  MOMENT,  INCH-POUNDS" ;M0 
70  J*0 

90  U*3. 14159265 
100  IF  A  <*  B  THEN  H*A/2 
110  IF  A  >  B  THEN  H*B/2 
120  Ql*<  H- . 5*TH)/B 

130  D*2»< ATN<  Q1/SQR< -Q1 *Q1  +  1  ) ) ) 

131  LPRINT  "  JAMMING  ANGLE,  RADIANS" ;D 

132  LPRINT  "  ANGLE;  MOMENT/MO;  MOMENT" 
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140  X=SU*<THA2>/4 


150  FOR  K=D/10  TO  H#D/10  STEP  D/10 

151  PI  =  K/2-KA 3/48+ KA 5/3840 

152  P2  *  1 -K A 2/8+ KA 4/384 

153  P3  -  B*P1 

154  P4  =  <2*H  -  P3> 

155  P5  =  B*P2 

156  P6  —  H*P1 

157  P7  =  P1/P2 

165  YB  -  P5-<  SQR<  P3#P4) ) 

170  R=  < . 0  7-K/70 ) *H 

180  ZA=<  P3*P1 ) -P6+  < SQR<  P3#P4> ) *P2 

190  Q1=1-<B/H)*P1 

200  Q2=ATN<Q1/SQR<-Q1*Q1+1 ) ) 

210  W1=2*A*<  U/2-K/2-Q2) 

220  UJ2=A*<U-2*Q2) 

230  W3=<P3*P1-P6+<SQR<P3*P4) ) #P2> *< U-2*Q2> 
240  U4=H*U 

250  YC=<H*P7+P5-<SGR<P3*P4> ) )/< 1+P7A2) 

260  XC=YC*P7 

270  Q1=ZA/<SQR<  <  <H-XC)A2>  +  <  <YC-YB> A2> > ) 

280  Q3*ATN<Q1> 

290  W5*2*B*Q3 

300  W7»4/3*<ZA/R)*<SQR<  < HA2> + < YBA2> + < ZAA2> ) ) 
310  U6»2*<H/R>*ZA 
320  Ql-ZA/YB 
330  Q3*ATN<  Q1 ) 
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340  W8=2*<A*<K/2)+2*H*Q3> 

350  J=J+1 

340  W< J>=UI1+W2+W3+W4+W5+W6+W7+W8 

370  NEXT  K 

380  FOR  l  TO  J-l 

390  Y=<TH/<B*<4*A+B>  > )*<W< IX+ 1 > -UK  IX> )/< D/1  0 ) 

391  M1=Y*M0 

400  LPRINT  ■  ■ ;IX*D/10jYjMl 

401  NEXT  IX 

402  LPRINT  ■  ■ 

403  LPRINT  "  " 

404  LPRINT  •  * 

410  LPRINT  “  * 

420  END 
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Appendix  C 


■  \  10  INPUT  -TEMPERATURE  OF  EXPOSURE,  DEGREES  F“ ;T 

30  INPUT  "FLANGE  WIDTH,  INCHES" ;A 
40  INPUT  "WEB  HEIGHT,  INCHES" ;B 
2  50  INPUT  "THICKNESS,  INCHES" ;TH 

61  LPRINT  "  " 

62  LPRINT  "  " 

•  70  LPRINT  "  TEMPERATURE  "  ;T 

30  LPRINT  "  FLANGE  WIDTH"; A 

90  LPRINT  "  WEB  HEIGHT" ;B 

S  100  LPRINT  “  THICKNESS" ;TH 

120  Tl=<A*B*TH+B*B»TH/4) 

130  T2=-20 . 5*T+400Q0 

135  IF  T>=1 951  THEN  T2=4.5 

137  M0=T2*T1 

138  LPRINT  "  FULLY  PLASTIC  MOMENT,  INCH-POUNDS" ;M0 

*  140  T3=0 . 36*3 .141 593*TH*( B+2*A> 

150  T4=l/<  B*<  B+4*A) ) 

170  LPRINT"  ANGLE;  MOMENT/MO;  MOMENT* 

B  171  IF  A<=B  THEN  J=A/2 

172  IF  A>B  THEN  J=B/2 
175  Ql=< J-.5«TH)/B 

*■  176  D=2*< ATN< Q1/SQR< -Q1 *Q1 ♦ 1 ) > > 

180  FOR  K=D/10  TO  D  STEP  D/10 


181  T5=TAN<K/2> 
190  M=T3*T4/T5 
195  Ml =M*M0 

200  LPRINT  ■  " ; K 

201  NEXT  K 

202  LPRINT  "  " 

203  LPRINT  "  * 

204  LPRINT  ■  " 


210  LPRINT 


n  n 
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